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Abstract 


We  examine  the  limited  aperture  problem  of  physical  optics 
farfield  inverse  scattering,  i.e.,  the  problem  of  identifying  a 
target  from  an  analysis  of  band-limited  viewing  apertures.  We 
show  (given  infcmation  from  all  directions)  that  the  target  may 
be  completely  identified  from  an  analysis  of  high  frequency,  band- 
limited  data.  If  the  directions  of  viewing  angles  are  limited  as 
well,  it  is  shown  that  the  target  surface  can  be  identified  where 
the  target  normal  lies  within  the  range  of  viewing  directions; 
outside  of  this  range,  no  information  is  available.  It  is  shown 
that  this  phenomenon  is  totally  a feature  of  the  Fourier  transform 
of  characteristic  (non-zero)  functions  and  independent  of  the 
inverse  scattering  formalism.  Numerical  examples  are  given  for  the 
case  of  a perfectly  reflecting  circular  cylinder 


1. 


Introduction. 


Physical  optics  inverse  scattering  constitutes  a technique 
whereby  the  shape  and  location  of  an  acoustical  or  electromagnetic 
scattering  target  may  be  determined  from  an  analysis  of  back- 
scattered,  far-field  data.  The  fundamental  identity  of  the  method 
(first  discovered  by  Bojarski,  [l]  ) relates  a quantity  p+(k) 
(that  is  proportional  to  the  back-scattering  cross-section  of 
radar  analysis),  to  the  Fourier  transform  of  the  characteristic 
function  of  the  scattering  target;  i.e.,  a function  which  is  unity 
inside  the  target,  and  zero  outside. 

The  key  assumptions  leading  to  the  Bojarski  identity  are 
(1)  that  the  target  may  be  modelled  as  an  acoustically  hard  or 
soft  reflector  (or  as  a perfect  conductor  in  the  electromagnetic 
case) ; (2)  that  the  physical  optics  approximation  is  valid  for 
the  back-scattered  fields;  and  (3)  that  the  source  and  observation 
ranges  are  large  compared  to  a typical  target  dimension.  The 
fundamental  identity*  which  is  derived  briefly  in  Appendix  A 
for  the  reduced  wave  equation,  states  that 

~ P+(jc)  + P*(-k) 

P(  k)  = = , (1.1) 

K2 

A* 

where  T(X)  is  the  target  characteristic  function,  and  r(£)  is 
its  Fourier  transform.  In  this  expression,  the  transform  vector 
jc  is  given  by 


k * 2kX  k = | , (1.2) 


Below,  we  shall  use  the  alliteration  POFFIS  for  "Physical  Optics 
Farfield  Inverse  Scattering". 
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where  k is  the  wave  propagation  number,  a)  is  the  (time-harmonic) 

A 

frequency,  and  c is  the  propagation  speed;  the  vector  X is  a 
unit  vector  pointing  in  the  direction  of  the  source-observation 
position.  The  symbol  (*)  denotes  complex  conjugation. 

A » 

Therefore,  if  r(ic)  could  be  measured  for  a sufficiently  large 
class  of  vectors  k_s  one  could  in  principal  determine  r(X)  by 
direct  Fourier  inversion.  As  (1.2)clearly  implies,  this  would 
necessitate  that  measurements  be  performed  at  all  frequencies  and 
at  all  aspect  angles  with  respect  to  the  target  about  which  informa- 
tion is  desired. 

An  immediate  difficulty  which  arises  with  the  method,  there- 
fore, is  the  lack  of  complete  information  in  transform  space;  the 
physical  optics  approximation  is  valid  only  in  the  high-frequency 
regime,  i.e.,  as  the  wavenumber  k approaches  infinity.  Further- 
more, data  that  would  be  available  from  an  actual  experiment  would 
most  likely  be  confined  to  some  band-limited  spectrum.  A further 
limitation  cf  the  POFFIS  method  occurs  because  data  can 
seldom  be  measured  at  all  aspect-angles  with  respect  to 
the  unknown  target. 

These  difficulties  lead  to  the  limited-aperture 
problem  of  the  FOFFIS  method  — the  problem  of  identifying 

A 

the  target  structure  from  the  limited  knowledge  of  T (k ) in  K-space. 
The  objective  of  this  paper  ■'s  to  present  solutions  to  the  limited- 
aperture  problem  for  a physically  reasonable  class  of  frequency 
and  aspect-angles  limited  apertures. 
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We  show  (given  information  from  all  directions)  that  the  target 
may  be  completely  identified  from  an  analysis  of  high-frequency, 
band-limited  data.  If  Lne  aperture  is  limited  in  aspect  angles  as 
well,  it  will  be  shown  that  the  target  surface  can  be  unambiguously 
determined  at  those  points  where  the  outward  normal  to  the  target 
surface  lies  within  the  range  of  directions  defined  by  the  aperture 
in  j<-space.  Outside  of  this  range,  no  information  is  possible  by 
this  method. 

The  limited-aperture  problem  has  been  previously  examined 
(from  an  analytical  standpoint)  by  Lewis  [2]  and  Perry  [3]  . Lewis 
considered  a number  of  aspect-angles  limited  apertures,  with  the 
(physically  unreasonable)  assumption  that  information  was  available 
for  all  frequencies.  Perry  concluded  that  the  absence  of  low- 
frequency  information  leads  to  an  ill-posed  problem  which  is  numer- 
ically unstable.  The  primary  difficulty  is  that  the  low-frequency 
data  constitutes  a principal  contribution  to  the  gross  resolution  of 
the  function  T(X)  . 

Additional  tests  of  Bojarski's  results  were  carried  out  numeri- 
cally in  papers  by  Tabbara  [4,  5]  and  Rosenbaum-Raz  [6],  Their  re- 
sults confirm  the  difficulties  with  resolution  due  to  bandlimiting. 
One  purpose  of  this  paper  is  to  show  how  this  shortcoming  of  the 
method  can  be  overcome. 

The  key  feature  of  the  method  we  describe  here  is  that  we  re- 
place the  problem  of  determining  T(X)  by  the  problem  of  determining 

A 

the  directional  derivative  of  T (X)  in  some  direction  £ • ^he 
directional  derivative  £ • Vr(X)  is  a highly  singular  function 
whose  essential  features  can  be  resolved  from  an  analysis  of  high- 

A 

frequency  data.  The  process  of  differentiating  in  the  £,  direction 


4 


in  X-space  corresponds  to  multiplication  by  £ • £ in  £-space;  in 
this  way,  one  simultaneously  suppresses  low-frequency  data  while  en- 
hancing high-frequency  data*  Indeed,  this  approach  was  first  sug- 
gested in  Bojarski's  initial  work  [1],  and  the  present  work  supports 
and  extends  his  ideas. 

In  section  (2)  are  displayed  a number  of  computer-generated 
examples  for  a variety  of  limited-data  apertures.  These  examples 
were  generated  from  the  exact  back-scattered,  far-field  data  which 
was  numerically  computed  according  to  a scheme  by  Greenspan  and 
Werner  [7]. 

We  remark  that  the  advantage  gained  by  calculating  the  Fourier 
transform  of  a directional  derivative  is  wholly  a feature  of  Fourier 
analysis.  It  has  nothing  to  do  with  Bojarski's  physical  optics 
identity,  nor  has  it  anything  to  do  with  the  discretization  and 
specific  calculation  technique  (FFT  in  this  work)  used  in  the  in- 
version of  the  transform.  It  is  shown  in  Section  3 by  asymptotic 
methods  why  the  band-limited  transform  behaves  as  it  does. 

More  precisely,  in  two  dimensions,  we  express  the  band-limited 
inverse  Fourier  transform  of  the  directional  derivative  as  a three- 
fold integral  in  arc  length  over  the  boundary  of  the  target,  and 
K and  0 , the  magnitude  and  angle  of  the  transform  vector  £ . 

Under  the  assumption  of  high  frequency  band-limiting,  the  first  and 
third  of  these  integrals  can  be  calculated  asymptotically  by  the 
method  of  two-dimensional  stationary  phase  [8].  The  integration 
in  magnitude  £ is  carried  out  exactly.  The  result  of  this  cal- 
culation explains  the  qualitative  features  of  the  numerical  results 

■k 

Multiplication  by  higher  powers  of  ic,  correspond  to  higher  derivatives, 
further  subdues  the  effect  of  low  frequence  data  and  further  enhances 
the  effect  of  high  frequency  data. 
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shown  in  Section  (2).  Furthermore,  stationary  phase  analysis  brings 
out  the  features  of  the  aspect-angles  problem,  as  well;  namely,  given 
a certain  range  of  aspect  angles  in  transform  space,  the  inverse 
transform  is  asymptotically  of  lower  order  outside  of  the  analogous 
aperture  in  physical  space. 

The  extension  of  this  analysis  to  the  case  of  three  dimensions 


is  discussed  in  Appendix  B. 
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2.  Numerical  Results 

Figures  (2.1-8)  display  the  results  of  the  method.  They  were 
generated  in  the  following  way.  The  reduced  wave  equation  was 
solved  for  a point  source  in  the  exterior  of  a totally  reflecting 
circular  cylinder  by  a method  due  to  Greenspan  and  Werner  [7].  The 
phase  and  range  normalized  backscattered  amplitude  (defined  below 
(A. 11)  in  the  Appendix)  is  calculated  to  five  place  accuracy.  This 
is  then  substituted  into  Bojarski's  basic  identity  (1.1)  to  yield 

A 

T(k)  . Table  I compares  the  exact  Fourier  transform  for  the  target 
cross-section  with  that  predicted  by  the  inverse  diffraction  tech- 
nique. Agreement  is  poor  for  wave  numbers  below  ka  = 10,  emphaiz- 
ing  the  need  for  high  frequency,  band-limited  analysis.  We  remark 
that  the  source-observation  range  was  25  radii,  indeed  in  the  far 
field,  so  that  this  could  not  have  accounted  for  the  discrepancy  at 
moderate  values  of  ka. 

The  backscattered  data  was  used  in  the  formula 

l * V r(X)  = (2tt)~2  j i l * £ TOO  e1— dK  , (2.1) 

for  the  directional  derivative  of  the  characteristic  function. 

Here  the  domain  of  integration  is  all  of  K-space.  In  particular, 
we  chose 

l = (1,0)  , (2.2) 

corresponding  to  transverse  differentiation  of  F(X)  , an  proceed  to 
calculate  (2.1)  over  finite  domains  representing  band-limited  and 
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aspected-angles  limited  apertures. 

The  inverse  transform  is  calculated  via  a fast-Fourier  trans- 
form (FFT)  routine.  Band-limited  data  around  k = 10,  20  and 
was  used,  with  a mesh  size  of  128  x 128  (of  which  most  entries  were 
zero)  and  mesh-widths  Ak^  = Ak^  =1-1  • The  resulting  function 
was  then  computer-displayed  to  yield  the  graphical  results  shown 
in  Figures  (2.1-7) . 

The  resolution  in  those  diagrams  is  seen  to  be  approximately 
one-half  the  minimum  wavelength  associated  with  the  domain  of  inte- 
gration in  K . The  clarity  of  the  results  is  also  seen  to  be 
directly  a function  of  the  frequency  band-width,  as  expected. 
Furthermore,  we  believe  that  better  resolution  may  be  attained  by 
curve  fitting  the  actual  output  to  the  results  predicted  by  the 
asymptotic  analysis  in  the  next  section.  The  "smearing"  effect 
near  the  lower  and  upper  portions  of  the  target  is  due  to  the  fact 
that  we  calculate  the  transverse  derivative  of  T(X),  and  this 
derivative  approaches  zero  at  these  limits.  This  minor  shortcoming 
maybe  overcome  in  one  of  the  following  ways:  (1)  calculating  a direc- 
tional derivative  which  does  not  vanish  in  these  ximits  (different 

A 

£ ) or  (2)  by  calculating  a higher  derivative  (for  example,  by 
2 A 

multiplying  by  k rather  than  £ * £ before  Inverting  the  Fourier 
transform) . 

Figures  (3.5-7)  show  the  results  from  a viewing  aperture  which 
is  limited  in  direction,  as  well.  The  theory  developed  in  the  next 
section  will  predict  good  resolution  of  the  target  surface  where 
the  outward  normal  to  the  target  lies  within  the  family  of  directions 


specified  by  the  viewing  aperture;  outside  of  this  range  no  infor- 
mation is  available  . 

One  further  remark  is  due  concerning  the  computer  plotting 
results.  The  theory  developed  in  the  next  section  is  based  upon 
leading  order  asymptotic  effects.  There  are  also  secondary  effects, 
corresponding  to  boundary  integral  contributions,  etc.,  but  these 
do  not  appear  due  to  normalization  in  the  plotting  process. 


"V 
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3.  Limited  Aperture  Analysis 


We  proceed  in  this  section  with  an  analysis  of  the  integral 
(2.1).  Let  us  denote  the  region  of  £-space  where  the  measured 

A 

function  T(k)  is  known  by  D , and  the  characteristic  function  of 
D by  A(]<)  : 


£<j0  = { 


1 

o 


< ED 
K $ D 


We  consider  the  function 


?X(K)  = i £ • K r(K)  . 


(3.1) 


This  is  the  Fourier  transform  of  the  directional  derivative  of 

A 

T(X)  in  the  £ direction.  We  note  that  the  product  of  r^0<) 

A 

and  the  aperture  function  A(jc)  , 

H(K,  £)  = AOc)  (ic)  , (3.2) 

is  known  everywhere  in  £-space.  We  consider  first  the  two-dimen- 

sional version  of  (3.2)  (Applicable  when  the  target  is  known  to  be 
a cylinder  whose  cross-sectional  profile  is  sought).  We  introduce 
as  new  coordinates  arc  length  a along  the  boundary  curve  3B 
and  distance  normal  to  the  boundary,  denoted  by  s (see  Figure 
(3.1)).  The  directional  derivative,  r^(X)  » may  be  represented  as 

rx(X)  - l • 7 r(X)  = l • n 6(s)  . 

Here  n = n(a)  is  the  outward  unit  normal  to  3B  . A Fourier  in- 

version of  (3.1)  and  use  of  the  convolution  theorem  for  Fourier 


(swjMrrr*"  — , 
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} 


transforms  leads  to 

H (X,£)  - (2ir)-2  n (a)  * £ A(X  - X*  (a))da  (3.3) 

3B 

The  objective  of  this  analysis  is  to  evaluate  0*3)  asymptoti- 
cally in  the  high  frequency  limit,  i.e.,  for  k = target  anywhere 
in  V . We  use  the  explicit  representation  of  A(X ) , 

A(X)  = e1- - d£  , 

D 

to  rewrite  (3.3)  as 

H(k,£)  = (2tt)~2  f da  n(a)  • £ ff  e1-  (-"-(a))dK  (3.4) 

J . 

3B  D 

The  three-dimensional  analogue  of  (3.4)  is 

H(X,£)  * (2tt)~3  | da  n(a)  • £ | d1-  (--X’(o))dic  , (3.5) 

3B  D 

where  a = (^>02)  represents  surface  coordinates  over  3B  (Figure 
(3.2)).  Analysis  of  (3.5)  is  carried  out  in  Appendix  B and  we 
proceed  with  the  analysis  of  (3.4)  here. 


Consider  an  aperture  which  is  limited  in  both  frequency  and 
directions,  as  illustrated  in  Figure  (3.3);  note  that  the  £-space 
aperture  D must  be  symmetric  about  the  origin  by  virtue  of  the 
analysis  leading  to  the  physical  optics  identity  (Appendix  A, 
equation  (A. 13) )• (When  the  aperture  directions  are  restricted  to 
one-sided  viewing,  it  is  possible  to  extend  the  data  symmetrically 
through  the  origin  in  £-space;  the  reconstructed  image  will  then 
be  symmetric  in  physical  space  (Lewis  [2]  ) , and  one  has  no  informa- 
tion concerning  the  opposite  side  of  the  target.) 

A A 

Define  j<  = kk  , where  k and  k are  the  magnitude  and 
directions  of  £ , and  rewrite  (3.4)  as 


r 
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(3.6) 

K 

H(x,o  y k<jk  y doMo) . t y .*;»>  • «-*•«'»» , 

K 3B  QD 

0 

the  angular  integration  in  this  expression  is  described  by  the 

a 

parameter  0,  which  is  taken  to  be  the  angle  of  the  vector  k. 

With  the  assumption  that  the  expression  in  (3.6) 

may  be  examined  via  a stationary  phase  analysis  in  the  variables 
(a,0).  With  the  phase  function  <J>(o,0)  defined  by 

$(0,0)  = k (0)  • (X  - X’(o))  , 

the  stationary  phase  conditions  become 

dX(o) 

<t>  ■ -k(0)  • T(o)  = 0 , T = , (3.7) 

0 da 

die(0) 

<t>  = x(0)  * (X-X’(a))  = 0 , t(0)  . (3.8) 

© da 

Let  E represent  the  portion  of  the  surface  3B  where  the  family 

A 

of  normals  n(a)  lie  within  the  range  of  directions  specified  by 

the  aperture  D.  Then  the  conditions  (3.7)  and  (3.8)  will  be  satisfied 

in  E when  the  vector  X - X' (a)  is  orthogonal  to  the  tangent 

A 

vector  T to  DB,  and  when  k(0)  is  collinear  with  X - X* (a) , 
as  illustrated  in  figure  (3.4). 

Since  the  jc_-space  aperture  is  symmetric  about  the  origin,  there 
exist  two  solutions,  0 = 0+,  for  which 

K(0+)  = -k  (0_)  . (3.9) 

There  may  also  exist  a multiple  of  solutions  in  the  arc-length 
variable  a;  however,  it  will  be  evident  from  further  analysis 


I 
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Chat  the  stationary  point  o = c for  which  (X  - X'(o))  is  a 

m — — 

minimum  will  dominate. 

Define  the  two-dimensional  Hessian  matrix  H by 

e J 


Hg  = det[$(cr,0) ] ; 


we  note  that 


dX'(cr  ) 2 

V -|x -x'(O|xc0  -| 2-1 


do 


In  this  expression,  y = X(°  ) is  the  curvature  of  the  target 

m 


surface  9B  at  o = a , and 

m 


!+l  X 

-i  X 


£ B 
{ B 


With  these  identifications,  H(X,5)  is  asymptotically  equal  to 


H(X,C)  * 


(9/ 


dxn(o 


_ ■ Cosf|x-X’ton)|} 
nr  5 . /TT~r7\  * 


yj |det(H  ) 


that  is, 


"<»■>  ' « (Sin(K,ix  - X'M>  - Sin^X  - X'M>) 

H(X>?)  * — 2==  < : : > (3.10) 


■rr\J|det(H ) 


X - X* (a  ) | 


m 


The  function  H(X, £)  thus  behaves  in  the  region  E (and  in  the 
high-frequency  limit)  as  a SinX/X  function  with  central  lobe  located 

A 

on  the  target  surface;  a trace  of  H(k, £)  would  therefore  yield  in- 
formation about  the  unknown  target  in  the  region  E. 

This  analysis  is  not  valid  when  the  Hessian  matrix  H becomes 

e 

singular;  the  vanishing  of  det(Hg)  occurs  on  the  evolute  of  the  target 
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surface  3B.  These  are  stationary  points  of  higher  order.  However, 
one  can  verify  that  the  amplitude  at  those  points  is  zero  and  hence 
the  contribution  from  such  a stationary  point  is  lower  order  than 
the  contributions  derived  above. 


The  calculations  leading  to  (3.10)  have  been  concerned  with 

an  estimate  of  contributions  arising  from  stationary  points  interior 

to  the  boundary  of  the  a - 0 domain  of  integration,  labeled  as 

3 (D,B)  in  Figure  (3.5).  (Note  that  D = S.  + S„+  . . . + S . 

i 2 8 

In  general,  contributions  from  stationary  points  on  the 
boundary  surface  3(D,B)  will  also  occur.  Define  the  vector  U. 
by 


U = 


V<$ 


|V<$|2  , 


V = 


\3o  ’ 30/ 


(3.11) 


and  note  that 


ge 


iK<E>  ie 


iw$ 


. U - i V • (UeiK^)  . 


— K 


Then  to  leading  order, 


H(X,§)  -v 


-iff/ 2 , 


(2ir)‘ 


/ dK  / 

kq  3(D,B) 


(3.12) 

(n(a)  • b( N • U)WeiK‘(0)  * ^^'^da 


where  a describes  the  one-parameter  surface  3(B,fi^), 


0 = 0(a)  , a = o(a)  , 


and 


The  unit  vector  N = N(a)  points  to  the  exterior  of  9(D,B). 

A stationary  phase  analysis  may  now  be  performed  on  the  remain- 
ing a-integral  in  (3.12);  these  points  of  stationarity  will  yield 
the  3 (D, B)  boundary  contributions.  The  objective  in  considering 
these  lower  order  contributions  is  to  show  that  no  information  is 
forthcoming  concerning  those  points  in  E on  3B  whose  outward 
normals  lie  outside  of  the  set  of  directions  implied  by  the  jc-space 
aperture  D;  outside  of  E,  no  interior  stationary  points  exist 
ana  the  boundary  terms  play  the  dominant  role. 

We  begin  by  considering  the  contributions  from  boundaries 
S5  and  S7;  a similar  analysis  will  yield  corresponding  results 
for  boundaries  Sg  and  Sg.  Furthermore,  the  total  contribution 
from  Sj  - can  be  seen  by  similar  methods  to  vanish  to  leading 
asymptotic  order. 

On  the  boundary  S5,  the  following  identifications  can  be 

made: 


0 = <J)Q  a = o 

=>j'Chw=i 

N = (0,-1) 


Equation  (3.12)  then  assumes  the  form 


-iTT/2  1 T 

H(X,§)  -v  2 f dK  f (n  • 0(N  • U)eiK^^o) 

(2tt)2  J J 

K0  °i 

Identify  the  phase  function  <J>5  as 


(X  - X'(o)) 


do. 
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$5  = K(V  * & " 2L’(o))5 


then  the  condition  of  stationary  phase  becomes 


deb 


da 


= -kU  ) • T (a)  = 0 , T = 


dX(a) 


da 


( .13) 


Condition  (3.13)  has  at  least  one  solution  a = a5,  where 


k(<*>0)  J-l(a5)  • 


The  factor  N • ]J  then  has  the  explicit  form 


N • U = 


-t(V  • (X-X'(°5>) 
(t(0q)  ♦ (X  - X’  (°5))]2 


T = 


die  (9) 
d0 


Note  further  that 


d^ 

da2 


dT(a  ) 

— - +X(a_), 

da 


where  x(o,-)  again  represents  curvature  at  the  point  a = a_;  for 
a convex  body,  this  quantity  is  positive.  A one-dimensional  sta- 
tionary phase  calculation  therefore  gives 

k. 


-ln/2 

H (X,5)  ^ — 

5 (x(0)l/2(2*)3/2 


1 


/ 


(N  • U)e 


i|K<(0o)  • (X  - X'  (O5))  + 


where  N_  is  to  outward  normal  in  (a  - 0)  space  to  S_. 

5 5 

An  analysis  of  the  boundary  S (where  0 = u + <£0)  will  yield 
a similar  result,  with  a stationary  solution  in  a again  at  a = a ; 

5 
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the  only  distinction  in  this  case  is  that  now 


k(u  + = -k(4>^) 


and 


(Ng  • U) 


- -(N?*  U) 


0 = 0n 


0 = TT  + 0 


The  sum  of  the  leading  order  contributions  from  Sg  and  s7> 
therefore,  gives 


A 


(n(qj)  • bdy  £(0^) 
(X(  o5))li/2( 27r) 3/2 


SinjKK(0o)  • (X  - X' (o5))  + tt/4| 
— 


dK, 


and  by  integration  by  parts  (retaining  only  the  leading  order  term) 


<n(a)  • 5)<V  U(0J 

h (X,0  * 2 ^ L 

5+7  ~ (X(  a))1/2(2Tr>3/2 


Cos J 1CK ( 0q)  • (X-  X'(C^))  + */*.} 
K 1 ^ ( 0^)  • (X  - X'  ( Og)) 


(3.14) 


The  principal  point  to  be  concluded  from  an  examination  of  (3.14) 


is  that  (in  the  exterior  of  the  region  E as  identified  in  Figure 

A 

(2.2))  H(k_,5)  exhibits  to  leading  order  a peaking  behavior  with 

central  lobe  occurring  in  the  vicinity  of  the  line 


where 


k(4>J  • T(o)  «=  0 ; 


however,  H(X,£)  gives  no  information  concerning  the  target  surface 
in  this  region. 


3.8 


Appendix  A. 

In  this  section,  we  present  a brief  derivation  of  the 
POFFIS  identity  for  acoustical  wave  propagation;  the 

electromagnetic  (vector)  case  has  been  outlined  by  Lewis  [3"]. 

As  indicated  in  figure  (A.l),  we  illuminate  the  acous- 
tically hard  or  soft  scatterer  B from  a time-harmonic  point 
source  located  at  Xq.  The  incident  field  is  thus  given 

by 


lyx.k) 


eiklx  - Xpl 
4*|X  ~ X0I 


(A.l) 


and  the  total  acoustic  field  U is  the  sum  of  the  incident  U 

T I 

and  a back-scattered  field  U , which  is  also  measured  at  the 

S 

point  XQ  : 


UT  " UI  + US 


(A.  2) 


The  total  acoustic  field  U is  required  to  satisfy  the  (time- 
reduced)  acoustic  wave  equation  with  source  at  X = h : 


(V2  + k2)UT  = -6(X  - X0)  . 


(A. 3) 


In  addition,  we  require  that  U satisfy  an  outgoing  radiation 
condition  at  infinity,  and  one  of  the  two  boundary  conditions 


3Ut 

— = 0 X e B (+)  , (A. 4) 

8n 


or 
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U = 0 X e B (-)  • 
T “ 


(A.  5) 


The  first  (+)  boundary  condition  models  B as  an  acoustically 
hard  reflector;  the  second  (-)  condition  is  the  acoustically 
soft  case.  Using  (A. 1)  - (A. 5)  and  the  requirement  that  U 

S 

also  be  outgoing  at  infinity,  one  may  write  a Kirchoff  integral 
expression  for  : 


Wk>  =/{vr.w 


300^  - X') 

3n' 


- G(*o  - x')- 


3Us(X\k) 
3n' 


|d  S' 


, (A. 6) 


3B 


where  the  integration  is  performed  over  the  surface  (denoted  3B) 
of  the  scatterer  B . In  this  expression,  G(X)  is  the  free-space 
outgoing  Green's  function  given  by 


G(X)  = 


4n  |xj 


(A. 7) 


We  now  make  a crucial  assumption  concerning  the  integrand  of  (A. 6). 
Under  the  integration  sign  appear  the  unknown  surface  field  U and 


the  normal  derivative 


au. 


lates  U’s  and 
the  scatterer  by 


3n 


3u_ 


3n 


S 


. The  physical  optics  approximation  re- 


to  the  incident  field  U on  the  surface  of 

I 


U 'v  ± U 
S I 


3U 

3n 


„ au 

S T 

^ + — i. 

3n 


3U? 

v, 

3n 


U = — §-  = 0 


) 


XcL 


(A.  8) 


X £ D 
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where  the  (±)  notation  has  the  same  connotation  as  before.  We 
define  the  lit  side  of  the  body  (denoted  L)  and  the  dark  side 
(D)  by  constructing  tangential  rays  from  source  to  target,  as 
indicated  in  figure  (A. 2).  The  envelope  of  the  set  of  points 
where  the  rays  strike  the  body  will  in  general  define  a closed 
curve  £ separating  the  body  into  lit  and  dark  sides. 

Strictly  speaking,  (A. 8)  constitutes  a valid  identification 
only  for  the  case  of  plane-wave  incidence  against  a plane-wave 
reflector  and  in  all  other  cases  constitutes  an  approximation. 

But  the  approximation  is  increasingly  valid  for  small  wavelength 
(large  k)  propagation,  where  locally  the  indicent  wavefront  and 
the  target  surface  may  be  approximated  by  their  tangent  planes  [9]. 
The  separation  of  the  surface  8B  into  a lit  and  a dark  side 
also  is  strictly  valid  only  for  convex  scatterers. 

We  may  now  use  (A. 8)  to  rewrite  (A. 6)  as 


us(x0,k)  * ± 


dS'  . 


(A. 9) 


Into  this  result,  we  now  introduce  the  far-field  approximations. 

Under  the  assumption  that  |Xq|  )£>  |x' | (distance  to  the  source- 
receiver  point  is  large  compared  to  a typical  target  dimension)  , 
we  write  |XQ  - X' ( a 3 two  terms  in  a binomial  series: 

1*0  - X’l  * X„  - X0  * X'  ’ (A. 10) 

a 

where  XQ  = { X0 I > and  XQ  represents  a unit  vector  in  the  direction 
of  X.  . With  these  identifications,  we  have 
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u (XMO  * 


^ikcxo  - x0  • x’) 

AtiX„ 


G(X0  - X’)  * 


,ik(X0  - X0  • X') 


4nXr 


and  equation  (A. 9)  becomes 


Ug(X,k)  'v-  ± 


e^kXp 

(4nX0)2 


dS’  . 


We  may  now  perform  the  normal  differentiation  to  obtain 


U (X,k)  a,  + .(■21k)e !!. 

S (4ttX0)  2 


y -i2kX  ■ 

Xne  0 


*’dS’ 


(A. 11) 


Defining  the  phase  and  range  normalized  backscattered  field 

p+(x0*k>  by 

i2kX0 

V*o’k>  -p+(x0»k)  > 

S ^ (AttXq)-  + 


then  (A. 11)  becomes 


p+(XQ,k)  = +(2ik) 


" -2ikX  « 
Xne  0 


(A. 12) 


/ 

The  domain  of  integration  in  (A. 12)  is  only  over  the  lit  position 
L of  the  total  surface  3B  . With  the  objective  of  obtaining  a 
closed  surface  integral,  we  now  write  the  analogous  expression 
when  the  source-receiver  point  is  reflected  through  the  origin 
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[Figure  (A. 3)]  . An  analysis  similar  to  that  leading  to  (A.l?) 
gives 


P_(-Xo,k) 


= ±(2ik) 


//: 


X0e-2ikV 


dS ' 


(A. 13) 


Note  that  the  lit  and  dark  sides  have  been  reversed;  this  assump- 
tion is  valid  only  when  the  origin  of  coordinates  is  selected  to 
lie  within  reasonable  proximity  to  the  target  compared  to  the 
distance  to  the  source-receiver  point. 

From  (A, 12)  and  (A. 13),  we  now  conclude  that 


p+(XQ,k)  + p*(-X0,k)  = +(2ik) 


and  an  application  of  the  divergence  theorem  leads  to 


P+(£o>10  + P_(-Xo>k) 


(2k) : 


-Iff 


-2ikX  • X'  1/s 

e 0 — dV'  . t.A.14) 


B 


Now  introduce  the  characteristic  function  of  the  volume  region 
occupied  by  B : 


r(x)  = 


1 XcB 


0 X i B 


note  that  the  Fourier  transform  of 


r(x) 


is 
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Appendix  B. 

In  this  section,  we  outline  the  analogous  calculations  leading  to 
the  three-dimensional  version  of  (2.10).  Let  D again  represent  an 
aperture  which  is  both  band  and  aspect  angles  limited,  as  shown  in 
Figure  (B.l).  Parametricly , D is  represented  by 


* e [d>0,  * - 0J  U [it  + d>0,  2tt  <f>0] 

0 e [0O,  it  - ©J  U [it  + 0Q,  2i,  - 0j 


where  $ and  0 are  the  polar  and  azimuthal  angles  of  k, 
respectively:  k = (Cosd>Sin0,  Sin<#>in0,  Cos0) . A separation  of  the 

magnitude  and  angular  portions  of  the  jc- integration  in  (2.3)  leads 
to 


H(X,£)  = fr)7  K2dlC  / don (a)  • 

K-  3B 


Choosing 


iK x(d>,0)  • (X  - X'  (£) ) 


Sin0  d(£d0 


D 


(B.l) 


<i>”  ic(<^,0)  * (X  - X 1 (a)  ) , 


the  stationary  phase  conditions  become 


^?=  t • (X  - X'(a))  = 0 
80  0 - - - 


(B.2) 


• (X  - X' (0))  = 0 


(B.3) 


r 


£$  = Sin0  x 
Z4>  <P 
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K(d>,0) 

3o 

1 

o 

1! 

*T 

f ’ 

• 

(B.4) 

<(<*>,©) 

3a 

0 

dX’(a) 

• 0 

(B.5) 

Conditions  (B.l)  - (B.5)  imply  that  the  unit  vector  k lies  in 
the  tangent  plane  to  3B.  A point  of  stationarity  in  (c  ,o  ) then 
occurs  where  X - X.*  ^es  a^on8  the  normal  line  to  3B.  A 

lengthy  calculation  then  leads  to 


H(X,0  * - 


£ • n SinK|x  - X' 


K = K. 


v |det(ujx  - X' |F  + G) | 1/2 |X  - X' 


(B.6) 


where 


!+l  X e B 

-1  X { B 


and  F and  G represent  the  first  and  second  fundamental  co- 
efficients of  3B.  Again  the  target  surface  3B  may  be  identified 

by  the  peaking  of  the  SinX/X  function  when  k > k » 1;  this 

1 0 

gives  information  only  about  the  portions  E of  the  surface  3B 
that  have  normal  directions  coinciding  with  the  set  of  directions 
contained  in  the  K-space  aperture  D, 

A calculation  similar  to  that  leading  to  (2.14)  again  leads  to  the 
conclusion  that  no  information  concerning  3B  is  available  outside  of 
the  E region. 
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